On self-protecting singlets in cuprate superconductors 

J. RohleJ] 

Universitdt zu Koln, Ziilpicher Str. 77, D-50937 Koln, Germany 
(Dated: February 2, 2008) 

Abstract 

The basal area (Cu-Cu grid) of the cuprate superconductors not only tends to shrink on hole 
doping, as expected from single electron quantum chemistry, but exhibits also an electronically 
incompressible "hump" around optimum doping Uopt — 0.16. The hump collapses near critical 
doping Uopt ~ 0.19. We analyze the origin of the hump in terms of a classical liquid of interacting 
incompressible particles in a container with antiferromagnetic walls. Oxygen holes interacting 
with the wall form singlets, protect themselves against other holes by an incompressible "spin 
fence" , and thus interact also with the lattice. Occupation of the Cu02 lattice with holes must 
therefore follow a non-double-occupant constraint also for the oxygen cage enclosing the copper 
hole. Closest packing of self-protecting singlets is found to occur around critical doping; closest 
packing paired self-protecting singlets around optimum doping. These singlet-states are bosonic, 
but are not magnetic polarons. 
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INTRODUCTION 



The lattice parameters of the cuprate superconductors are well known strong functions 
of doping. It is however only poorly understood whether their doping dependence is ex- 
trinsically controlled from the complex stereochemistry of the doping blocks outside, or 
intrinsically by the properties of the strongly correlated electron states residing inside the 
Cu02 planes. The external chemistry and the embedded metallic layer reside in the same 
unit cell, and hence both are expected to co-determine the doping induced variations of the 
cell parameters. How to disentangle the most important intrinsic effects from the extrinsic 
ones? Detailed inspection of many of the available lattice data , displaying 

doping dependencies in the single- and multi-layer compounds of different families, shows 
forces of different nature to govern the behaviours of the a-, 6-parameters in the planar 
directions, and of the c- parameter in the perpendicular direction, a and b tend to contract 
with increasing hole concentration, almost independent on the number of CUO2 layers and 
the chemistry of the doping block. Uniaxial strain from the dopants may mask or even 
invert the shortening of one of the basal bondlengths, e.g. in underdoped YBa2Cu30a; where 
a contracts as b expands But the basal area B, defined by the square of the basal Cu-Cu 
distances, turns out almost unaffected by strain from the doping block. Thus B allows better 
for comparisons between compounds from different families than the individual basal bond 
lengths. The c-parameter however seems to behave arbitrarily: as the basal plane contracts, 
c may expand, contract, or both 

mm 

. Apparently the problem of entangled intrinsic 
and extrinsic doping effects on the lattice seems to solve quite naturally: the behaviour of 
the basal plane is governed by the nearly-2D quantum liquid residing in it, while that of the 
perpendicular parameter is dominated by the chemistry of the environment. In this article 
we will focus on the doping dependence of the basal plane B{x), in particular on the origin 
of the strong hump around optimum doping, see Figs. [H |2l El IH 



EXPERIMENTAL DATA 

Figs, m 121 El in display the basal areas B of typical one- and two-layer cuprates as a 
function of doping. B refers to the grid of nn planar copper atoms as reported from x-ray or 
neutron diffraction measurements at room-temperature by Radaelli et al. j|2| (La2-a;Sr2,Cu04, 
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Tcmax = 36 K), Fukuoka et al. [3] (HgBa2CuOx, T^max = 96 K; and HgBa2CaCu20x, Tcmax = 
127 K) , B6ttger et al. [i], and Kaldis ^1] (Yi.j^CayBazCugO^, T„ = 92 K). The thick 
drawn out hnes connecting the data points are guides to the eye. 

The overall behavior of B{x) exhibits surprisingly strong similarities in all systems under 
comparison: i. As expected from the increasing covalency with hole doping the basal areas 
shrink between the insulator-metal transition and the strongly overdoped regime. The 
contraction is in the order of 1.5% ii. B{x) exhibits a hump centered around optimum 
doping. The hump is weakest in La2_a;Sra;Cu04, and strongest in HgBa2CaCu20a;. Its 
maximum is centered aroud optimum doping Xopt- 

Following quantum chemical approximations hole doping will remove electrons from the 
antibonding a*Cu3d^2_y202p^ y band, increase the amount of covalent character in the Cu- 
O bonds, and will shorten them. The hump indicates a significant deviation from this one 
electron bandstructure picture. The thin dashed straight lines in Figs. [H El IH are fitted 
to the data points at the strongly under- and overdoped ends and, extrapolated towards 
the underdoped-overdoped phase boundary, turn out to intersect around Xopt- They may 
serve as coarse approximations for the quantum chemical "background" Bo{x). Its change 
of slope around Xopt might indicate that the Cu-0 bonding changes from semicovalency in 
the underdoped to covalency in the overdoped regime. 

Evidently the hump stems from the strong correlations of the holes in the Cu02 planes. 
Its location at optimum doping points also to a connection with the superconductivity oc- 
curing at much lower temperatures than 300 K. Consider the compressibility of the quantum 
liquid in the planes is approximated by oc —d{B — Bo)/dn, where n is the number of 
holes/Cu/unit area, and n oc x. Then — close Xopt, i-e the quantum liquid is incom- 
pressible around the optimum hole concentration Uopt = 0.15 — 0.16, notably not at quarter 
filling n = 0.25. 

In some compounds the collapse of the hump occurs together with subtle structural 
instabilities in the crystallographic cell, of martensitic type in YBa2Cu30a: 0, ^, or of 
order-disorder type in HgBa2Cu20a; 0|. Evidently these lattice instabilities are connected 
with the transition into the overdoped regime, but are most likely not at its origin. 

Stable long range ordered nano domains might be another source of the hump. 
La2-xSr2;Cu04, due to its octahedral tilts the system most susceptible to stable nano do- 
mains, exhibits however the weakest hump. On the other hand the Hg-cuprates with their 



3 



almost flat Cu02 lattices, devoid of structural compliance, exhibit the strongest humps. 
THE MODEL 

Consider B{x) reproduces the equation of state of the quantum liquid conflned to the 
Cu02 lattice. Suppose that the temperature is sufficiently high and the density of classical 
particles sufficiently low, then B{x) exhibits a striking similiarity with the van der Waals 
equation, or related types describing classical real gases. Here we describe the particles by 
singlets of oxygen hole and copper spins. We address the question, how a quantum liquid, 
created through singlet formation in a hole doped Cu02 lattice, may be connected with 
lattice properties such as B{x) and its hump around Xopt- Further down we will also try to 
justify why the problem may be phrased in terms of a classical Bose gas scheme. 

We consider the antiferromagnetic lattice of the copper spins as the walls of the "con- 
tainer" enclosing the available area for the spin singlets created upon hole doping. The 
moving holes exert a mean pressure (positive or negative) on the antiferromagnetic lattice 
given hj p oc nt where t is the average kinetic energy of the hole. Note that a Fermi liquid 
confined in the volume of its metal exerts its mean pressure only on the atomic lattice. 

A doped hole at an oxygen site destroys the antiferromagnetic order in its vicinity and 
thereby leads to an attractive interaction when it shares the region of depressed antifer- 
romagnetic order with a copper hole. The attractive interaction creating the spin singlet 
thus tends to keep the oxygen and Cu atoms closer together than it would be the case for 
noninteracting holes. Bottger and Dichtel P| find from a three-band Hubbard model an 
oxygen displacement per hole by -0.04 A along the Cu-0 bond. The creation of singlets has 
thus the same effect as a slight compression of the basal area complying with an increasing 
pressure. Thus the oxygen holes interact also strongly with the atomic lattice "eating" spins 
in the antiferromagnetic wall of its container. 

On the other hand short-range repulsive forces between the singlets keep them sufficiently 
apart to prevent them from occupying the same places at the same time. The area occupied 
by a singlet themself must be thus subtracted from the area available to any other singlet 
in the container. For sufficiently high doping the repulsive interaction between singlets will 
start to outweigh the contraction of lattice driven by singlet creation. Long range attractive 
interactions between the singlets will pack them more closely in the container and thus may 
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act as an additional pressure. 

The potential energy U of the interaction between the singlets may be expressed only by 
their relative separation R, in its simplest form as 



. X 1 cx) : R < Ro 

U{R) = { (1) 

-Uo{Ro/Ry : R>Ro 

Rq is the minimum possible separation between incompressible singlets. The exponent s 
is e.g. ~ 6 in typical van der Waals gases. 

Long range interactive forces may create also paired singlets excluding a larger area than 
two unpaired singlets. Then the area available to other singlets will be even more reduced 
than by closely packed unpaired singlets. As a result formation of paired singlets will 
render the Cu02 lattice even less compressible than most closely packed unpaired singlets. 
It is therefore suggesting to assume that the optimum doped incompressible Cu02 lattice 
accomodates its holes {uopt = 0.16) in closest packed paired singlets. 



The self-protecting singlet (SPS) 

How large will be the incompressible areas covered by a paired singlet, and a unpaired 
singlet, respectively ? Paired singlets in the most closely packed conformation have to 
match the maxium of the hump at Uopt — 0.16. Unpaired singlets in the most closely 
packed conformation are expected to match a critical hole concentration Ucrit, slightly larger 
than Uopt. For n > n^pt paired singlets will start to overlap with each other and thus will be 
broken, and the hump will collapse. We locate the collapse of the hump around ricru — 0.19. 
For n > ricrit even unpaired singlets will be broken, but may create a new type of singlet 
state excluding a smaller area than the singlets at n < ridt. 

We consider a perfect antiferromagnetic Cu02 lattice. The non-double-occupant con- 
straint for the Cu3d^ sites requires that doping must create additional holes on the oxygen 
sites. Hirsch showed that without flipping Cu spins a doped oxygen hole can only move 
within a "cage" of four oxygen atoms surrounding the nearest Cu atom. This is due to 
the pha se coherence in the symmetric combination of the four oxygen states. Zhang and 
Ricejll| worked out that the binding energy of the resulting spin-singlet state of the sym- 
metric oxygen hole and the Cu hole is 4 times higher compared to that of a spin-singlet state 
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of an oxygen hole sitting at a fixed site, and the Cu hole. Thus the so-called Zhang-Rice 
(ZR) singlet distinguishes itself from other possible singlet states by its extraordinarily high 
stability. The cage determines the area covered by a ZR singlet and is given by dQ_Q = 
Here do-o is the nn oxygen-oxygen distance, and a the Cu-Cu distance. The ZR singlet 
has no magnetic interactions with all other copper holes. But two neighbored holes trying 
to share their oxygen cages with each other will feel a strong repulsion. Thus formation of 
nn singlets in connected cages will be very unlikely. Rather the ZR singlets will strongly 
repel each other, and thus will create an excluded area around their cages. 

To protect the symmetric oxygen hole in the cage from other holes, the excluded area has 
to extend at least over the four neighbored cages containing the 4 nn Cu spins. The 4 nn Cu 
spins are exempted from the formation of own singlet states but act as a "fence" enclosing 
the ZR singlet. We label this as "self-protection" of the ZR singlet for a safe life, at least 
within its spin fence. Hence a CusOig cluster occupied by one oxygen hole in the central 
cage appears to be the minimum conformation for a stable singlet, that is a self protecting 
singlet (SPS). It covers nine nn oxygen-oxygen squares with an area 3dQ_Q = 9/2a^. Fig. 
displays a SPS in the antiferromagnetic Cu02 lattice. The circle enclosing the central cage 
indicates the ZR singlet, the hatched cages the excluded area with the spin fence. 

The CU5O16 cluster comprising a self protecting singlet (SPS) looks alike a polaron, but 
is not a conventional spin polaron polarizing the Cu lattice ferromagnetically. Here holes 
and spins are in the same band unlike the situation in the magnetic polarons. 

Upon propagation of the SPS from a given site to a neighbored site the oxygen hole 
has to flip its spin, and will thus be able to create a singlet with the nn Cu hole. The 
Cu spin of the abandoned site will recover and hence a possible "loop-hole" in the spin 
fence will be closed. Thus the spin fence moves together with the propagating hole while 
the antiferromagnetic lattice remains intact. This scenario has appealing similarities with 
the "spin-bag" mechanism proposed by Schrieffer et al. in that both, the spin-bag and 
the self protecting singlet, are both polaron-like, but are not conventional spin-polarons 
with heavy masses. We understand however that a fundamental difference occurs in that 
in the spin-bag mechanism the propagating hole and its surrounding bag act as fermionic 
quasiparticle, wheras a hole propagating in a self protecting singlet has to be considered as 
a bosonic particle. 

We may phrase the description of the SPS in terms of another non-double-occupant 
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constraint for the doped holes: the grid of the 5 oxygen corner hnked oxygen cages in a 
CU5O16 cluster may be only singly occupied. Or alternatively expressed: the spin fences of 
self-protecting singlets are not allowed to overlap. 

Respecting both non-double-occupant constraints, that for the Cu sites and that for 
the oxygen cages, it is apparent that the maximum number of singlet states in a doubly 
constrained CUO2 lattice will be much smaller than in the only singly constrained case. We 
will show further below that the doubly constrained doping of holes in a CUO2 lattice will 
be generically inhomogenous for sufficiently high concentrations. Ignoring the non-double- 
occupant constraint for the oxygen cage will however allow for homogenous hole distributions 
at optimum doping, as e.g. in the RVB models jisj ]. 

Using the information from the experimental B{x) we will estimate the minimal size of 
the particles occupying the CUO2 lattice in the different regimes of doping. 



The paired self-protecting singlet (PSPS) 

Fig. ini displays two self-protecing singlets sharing a common oxygen site along the direc- 
tion of the Cu-0 bond. This conformation may lead to an exchange of holes between two 
intact ZR singlets (closed circles) through opened spin fences as indicated by the dashed 
circles. The thus connected SPS are fully antisymmetric in respect to the connecting oxygen 
and may form in principle a paired ZR singlet state as proposed by Miiller Self- 
protection of the "paired" singlet will be achieved by an enlarged spin fence as indicated 
by the cross-hatched squares in Fig. IHl The area, covered the by the spin fence of the 
PSPS, is nominally 12a^, significantly larger than 9a^, the area protected by the spin fences 
of two nonbonding SPS. Thus formation of a PSPS will reduce the area available to other 
singlets much more than most closely packed SPS. Notably the PSPS extends over 4a ~ 15 
A (a chain of four corner-linked oxygen cages), nearly coinciding with the experimentally 
established planar superconducting coherence length C^ab- 

Because of the the contraction of the Cu-0 bonds upon singlet formation both, the SPS 
and the PSPS, are expected to exhibit important vibrational properties. Notably the oxygen 
holes are expected to breathe within the rigid 2a cell of the spin fence as observed in e.g. 
optimum doped Lai.85Sro.i5Cu04 by McQueeney et al. jisl ]. 
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Characteristic distributions of SPS and PSPS 

Consider the Cu02 lattice as an array of corner linked oxygen cages each enclosing one Cu 
atom (Figs. 171 11111) . One SPS covers five of these cages and thus forms a 3 x 3 supercell. The 
3x3 supercells will create a checkerboard lattice of black (dark/red) fields with "crosses" of 
5 cages, and white (grey) fields with "squares" of 4 cages. Doping will distribute the holes 
as SPS over this checkerboard lattice. At sufficiently high concentrations the SPS will touch 
each other, form clusters, or most likely PSPSs. A full occupation of the lattice may be real- 
ized by many different distributions leading to a remarkable sequence of characteristic hole 
concentrations. Figs. ITIITUI sketch four possible distributions relevant for the understanding 
of the phase diagram of the cuprates. The SPS are indicated by objects with a circular cage 
(the singlet hole) enclosed by a square (the spin fence). 

Stripy long range ordering: n = 0.11 ~ 

In Fig. [7|each "cross" (dark/red) is occupied with one SPS. Hence the "squares" (grey) 
will not be available for the SPS. The SPS will thus align in crossing chains along a. One 
hole rich chain alternates with two empty chains creating therewith a 3a superstructure. 
Presumably pairing of neighbored SPS to PSPSs will be suppressed by long range ordering. 
Only 0.11 holes/Cu are accomodated by this "stripy" distribution. 

Critical Doping: n = 0.20 ~ Ucrit 

To achieve closer packing of holes the SPS have to be distributed over the crosses and 
the squares of the checkerboard lattice. Fig. |H1 displays a possible distribution for a high 
density of holes. Cells excluded by the spin fences are cross-hatched. The distribution for 
the maximum number of SPS will be realized by most closely packing at Ucrit- For n > Udt 
the spin fences will overlap and destroy the self-protection of the singlets. Interestingly this 
leads generically to an inhomogeneous hole distribution, apparently because the "crosses" 
and "squares" in the checkerboard lattice are not equivalent for the cross-shaped SPS. The 
inhomogenous hole distribution favors the formation of PSPSs as it is indicated by the cigar- 
shaped spin fences enclosing the two neighbored SPS connected along a. The SPS may only 
disappear for dopings exceeding ricrit- Boundary effects destroying the self-protection of the 
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singlets will result from the finite size of the crystal. SPS damaged by boundary effects 
are indicated by filled circles with dashed spin fences. About ~ 0.18 holes/Cu may be 
accomodated in this distribution. Closest packing of SPS accomodates 1/5 = 0.2 holes/Cu. 

Optimum Doping; n ~ 0.16 ~ Uopt 

Fig. ini displays as an example a distribution of nearly closest packed paired self protected 
singlets. About 0.15 holes/Cu may be accomodated by this distribution. The boundary 
effects are even more serious for the PSPSs than for the SPS: PSPSs will be destroyed 
within 3a from the edges of the crystal. Closest packing of PSPSs accomodates 1/6 = 0.166 
holes/Cu. 

"Destructive" doping: Udes > 0.2 

For rifles = 0.22 each SPS must overlap its spin fence with that of the nn SPS (see Fig. 
inH) . At such high hole concentrations the singlets will be unable to protect against other 
holes, break apart or possibly localize. Phrased in terms of the classical theory of gases: for 
rides the antiferromagnetic lattice may no longer serve as the walls of the container enclosing 
the liquid. 

SUMMARY AND CONCLUDING REMARKS 

The ubiquitous observation of a hump in B{x) of cuprate superconductors has led us to an 
attempt analyzing the quantum liquid confined to doped Cu02 lattices in terms of the theory 
of classical real gases or liquids. Most importantly here the walls of the container enclosing 
the quantum liquid are formed by the antiferromagnetic lattice which couples to the atomic 
lattice by creation of self-protecting spin singlets. It is basically internal magnetostriction, 
although of a novel type, that allows for insights into the structure of strongly correlated 
electrons from an analysis of the lattice parameters. 

We have shown that the hump is strong evidence for the existence of relatively large 
incompressible particles interacting with each other. The area of the particles can be es- 
timated from the location of the hump around x^pt = 0.16, and from its collapse around 
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Xcrit > Xopt, respectively. At a microscopic level neither simple (nn) singlets, nor the ZR- 
singlet will match the area of 6.25a^ to be covered for hard core repulsion at Uopt = 0.16. 
The area of pairs of PSPS however is close to this requirement. On the other hand the 
area 9/2a^ covered by an unpaired SPS leads to close packing at ricru — 0.20, and seems 
to match even the number of destructive doping rides — 0.22 where in IR experiments the 
bosonic excitations disappear 

The SPSs are most likely bosonic, and not fermionic. Clearly this and other issues of 
the SPS and the PSPS will raise a bunch of challenging questions, which certainly cannot 
be answered from our so far almost only phenomenological approach which lead us via a 
geometrical construction to these particles. It is clear however that displacive lattice degrees 
of freedom must be involved in the formation of the singlet states. Otherwise the hump would 
not occur. 

I appreciate stimulating and enlightening discussions with K. A. Miiller, T. Egami and 
E. Stoll during this conference. 



Electronic address: labbl2@uni-koeln.det URL: http : //www . uni-koeln . de/~abbl2 

E. Kaldis, in Handbook on the Physics and Chemistry of Rare Earths, edited by K. A. Gschnei- 

dner Jr., L. Eyring, and M. Maple (North-Holland, Amsterdam, 2001), vol. 31. 

R G. Radaelli, D. G. H. A. W. Mitchell, B. A. Hunter, J. L. Wagner, B. Dabrowski, K. G. 

Vandervoort, H. K. Viswanathan, and J. D. Joergensen, Phys. Rev. B 49, 4163 (1994). 

A. Fukuoka, A. Tokiwa-Yamamoto, M. Itoh, R. Usami, S. Adachi, and K. Tanabe, Phys. Rev. 

B 55, 6612 (1997). 

G. Bottger, I. Mangelschots, E. Kaldis, P. Fischer, C. Kriiger, and F. Fauth, J. Phys.: Condens. 
Matter 8, 8889 (1996). 

Q. Huang, J. W. Lynn, Q. Xiong, and C. W. Chu, Phys. Rev. B 52, 462 (1995). 

J. D. Joergensen, O. Chmaissem, J. L. Wagner, W. R. Jensen, B. Dabrowski, D. G. Hinks, 

and J. F. Mitchell, Physica C 282-287, 97 (1997). 

E. Kaldis, J. Rohler, E. Liarokapis, N. Poulakis, K. Conder, and P. Loeffen, Phys. Rev. Lett. 
79, 4894 (1997), cond-mat/9707196. 

J. Rohler, in Concepts in Electron Correlation, edited by V. Zlatic and A. Hewson (Kluwer, 



10 



Dordrecht, 2003), vol. 110 of NATO Science Series II, p. 55, cond-mat/0210560. 
[9] T. Bottger and K. Dichtel (2003), private communication. 
[10] J. E. Hirsch, Phys. Rev. Lett. 59, 228 (1987). 
[11] F. C. Zhang and T. M. Rice, Phys. Rev. B 37, 3759 (1988). 
[12] J. R. Schrieffer, X.-G. Weng, and S.-C. Zhang, Phys. Rev. Lett. 60, 944 (1988). 
[13] P. W. Anderson, The Theory of Superconductivity in the High-Tc Cuprates, Princeton Series 

in Physics (Princeton University Press, Princeton N.J., 1997). 
[14] K. A. Miiller (2003), this conference. 

[15] R. J. McQueeney, Y. Petrov, T. Egami, M. Yethiraj, G. Shirane, and Y. Endoh, Phys. Rev. 

Lett. 82, 628 (1999). 
[16] T. Timusk (2003), this conference. 



11 



CM 

to 

03 



03 
w 

CO 

OQ 



14.4 



14.35 - 



14.3 - 



lo 14.25 - 



14.2 - 



14.15 



1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
La S r CuO 




2-x X 4 — 


~ \ 

\ 


295 K : 


I underdoped 
~ 


T overdoped I 


cmax 

1 , , , , 



0.05 0.1 0.1 5 0.2 0.25 0.3 0.35 

x(Sr) 



FIG. 1: Basal area of one-layer La2_xSr2:Cu04 as a function of doping using neutron diffraction 



data of Radaelli et aL l2l- 



a 



15.15 



^ 15.1 - 



CO 

CD 



CO 
CO 

CO 
GQ 



15.05 - 



14.95 



1 1 1 1 1 1 1 1 1 1 1 1 1 
- underdoped T 

cm 


1 1 1 1 1 1 1 1 1 1 
overdoped 

ax 


7 HgBa^CuO^ 




~- 300 K 




1 , , , 



4.02 4.04 4.06 



4.08 
X 



4.1 4.12 4.14 



FIG. 2: Basal area of one-layer HgBa2CuOa; as a function of doping using x-ray diffraction data of 
Fukuoka et al. 



12 



CO 



< 

w 
CO 



14.92 
14.9 
14.88 
14.86 I- 



0.1 



0.2 



14.8 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

- underdoped T overdoped 

cmax 


• • * 


- 


• • . VvJ • 










T • — 






= Y Ca Ba Cu 




— 1 - y y 2 3 X 


• ^^'^^ • — 


: 300 K 







• 



6.2 6.4 



6.6 
X 



6.8 



G. 3: Basal area of two-layer Yi_j^CayBa2Cu30x using x-ray diffraction data of Bottger et al. 



(overdoped by Ca, y ~ 0.96 = const. ), and Kaldis and collaborators Q (underdoped). 



14.94 




X 



FIG. 4: Basal area of two-layer HgBa2CaCu203; as a function of doping using x-ray diffraction 
data of Fukuoka et aL 31 



13 



FIG. 5: Self-protecting singlet (SPS) in an antiferromagnetic Cu02 lattice. The open circle in the 
center of the hatched supercell of 3 x 3 nn oxygen-oxygen cells indicates the symmetric oxygen 
cage accomodating a ZR singlet of the copper and the oxygen hole. The cross-hatched squares 
protect the ZR singlet against other holes with their intact spins ("spin fence"). The SPS covers 
an area of 9/2a^ in the Cu02 lattice. 



FIG. 6: Two self-protecting singlets sharing a common oxygen site along the Cu-0 bond direction 
create one paired self-protecing singlet (PSPS). The spin fence in the cross-hatched squares protects 
the exchange between the two SPS against other holes and the loss of symmetry. The PSPS covers 
an area of 6a^ in the Cu02 lattice. The PSPS extends over four oxygen cages along a. 
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FIG. 7: "Stripy" doping: Perfectly ordered distribution of oxygen holes over the antiferromagnetic 
Cu02 lattice by self-protecting singlets indicated by the open circle and the squared spin fence. 
Note the checkerboard lattice of "crosses" (dark/red) and "squares" (grey). Only 0.11 holes/Cu 
are accomodated. 
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FIG. 8: Critical doping: Possible high density packing of oxygen holes in the antiferromagnetic 
Cu02 lattice by SPS and PSPSs over the "crosses" and "squares". Cross-hatched cells indicate 
protecting areas. Note that the distribution must is inhomogeneous favoring the formation of 
paired singlets (large cigar-shaped spin fences). Boundary effects are indicated by filled circles and 
dotted broken spin fences. In this example about 0.18 holes/Cu are accomodated. 
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FIG. 9: Optimum doping: Nearly closest packing of oxygen holes in the antiferromagnetic CuO 
lattice by PSPSs. Note that closest packing of PSPSs will pair all available SPS. In this exampl 
about 0.15 holes/Cu are accomodated. Boundary effects are indicated by dashed spin fences. 
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FIG. 10: Destructive doping: In this example 0.22 holes/Cu are accomodated. All spin fences are 
broken. 
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